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Abstract 



This works extends the recent study on the dielectric permittivity of crystals within the 
Hartree model [E. Cances and M. Lewin, Arch. Rational Mech. Anal, 197: 139-177, 2010] to 
the time-dependent setting. In particular, we prove the existence and uniqueness of the nonlin- 
ear Hartree dynamics (also called the random phase approximation in the physics literature), 
. in a suitable functional space allowing to describe a local defect embedded in a perfect crys- 

tal. We also give a rigorous mathematical definition of the microscopic frequency-dependent 
(— I , polarization matrix, and derive the macroscopic Maxwell-Gauss equation for insulating and 

" ■ semiconducting crystals, from a first order approximation of the nonlinear Hartree model, by 

means of homogenization arguments. 



I 1 Introduction 

\^ , A material subjected to a time-dependent perturbation usually does not respond instantaneously. 

Consistently with the causality principle, the linear response of the material can be expressed 
as the time convolution of some causal response function with the applied perturbation. The 
response properties are therefore frequency-dependent in general. This is the case for instance 
for the dielectric permittivity of the material, which allows to describe the linear response of the 
electronic density in terms of an applied external electric field [1] [23]. For molecules, a dipole 
moment is created, while for solids a more global charge redistribution, with possibly screening 
effects, occurs. 

For molecules, a convenient model to approximate the many-body quantum dynamics of the 
■ system is the time-dependent Hartree- Fock model, whose well-posedness is studied in |11[71[H]- In 

^ , the density matrix formulation of the Hartree- Fock model considered in [7] , the state of the system 

at time t is described by a density matrix 

7(t) e SiL^W')), < 7(0 < 1, (1) 

where S{L'^{S?)) denotes the space of bounded self-adjoint operators on L^(R'^), and where, 
for A and B in S{L'^{R^)), A < B means < [i) , B%l^) ^2 for aU -0 e L'^{m.^). To 

simplify the notation, we omit throughout this article the spin variable. This does not modify 
the mathematical structure of the equations. The condition < 7(t) < 1 is a translation of 
the Pauli exclusion principle in the language of one-body density matrices: two electrons cannot 
be at the same time in the same quantum state. The density matrix ^{t) is in fact trace-class 
since there is a finite number of electrons in the system (recall that in non-relativistic quantum 
mechanics, the trace of '){t) is stationary, and equal to the number of electrons). Therefore, an 
electronic density p^i^t) G L^(M'^) can be associated with the operator ^{t) [H]. We consider here 
the dynamics of the finite system within the time-dependent Hartree (also called time-dependent 
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reduced Hartree-Fock) approximation: 



and the dynamics is modified as 



with 

H° = --A + Kuc+fc(P7). 

where Vnuc is the potential generated by the nuclei, and 

is the Coulomb potential generated by the charge density p. The time-dependent Hartrce model 
is obtained from the time-dependent Hartree-Fock model by discarding the exchange term. It 
can also be viewed as the simplest model derived from time-dependent density functional theory 
(TDDFT, see for instance [15]), corresponding to the case when exchange-correlation is neglected. 
When an external perturbative potential v{t) is considered, the Hartree Hamiltonian reads 

H:; = -^A + Vnuc + Vcipj) + V, 

i^ = [i^.",7]. (2) 

The well-posedness of such dynamics is studied in [2] . Recently, extensions of the time-dependent 
Hartree-Fock models have been studied, in particular time-dependent multi-configuration mod- 
els lain]. 

Crystals arc infinite periodic assemblies of nuclei surrounded by their electronic clouds. The 
currently most popular models to approximate the dynamics of their electronic structures rely 
on TDDFT, and read as self-consistent nonlinear mean-field models. However, linear empirical 
models are sometimes used. In both linear empirical models and self-consistent nonlinear mean- 
field models, the electronic state of the crystal at time t is described by a one-body density matrix 
7(t) still satisfying ([1]). On the other hand, since there are infinitely many electrons in a crystal, 
7(t) is not trace-class. 

In linear empirical models, the electrons in the crystal experience an effective potential and do 
not interact with each other (except through the Pauli principle). In such models, a perfect crystal 
with periodic lattice TZ is characterized by a periodic Schrodinger operator H^^^. = — -f Vpev 
where Vpor is a given 7^-periodic effective potential. Provided Vpcr S Lf^^{M.^), i/p^r defines 
a self-adjoint operator on L^(R'^) with domain _ff^(M^), bounded from below, with well-known 
mathematical properties. In particular, the spectrum of -ffp^,. is purely absolutely continuous and 
composed of a countable number of (possibly overlapping) bands [22]. The ground state of the 
system is described by the one-body density matrix 

/per -^( — oo.epl V-'-'pcr/i 

where the real number ep, called the Fermi level, controls the number of electrons per unit cell. 
Here and in the sequel, 1b denotes the characteristic function of the Borcl set B C M.. Loosely 
speaking, the electrons fill the energy levels of -ffp^r to £f, and filling the N lowest energy bands 
amounts to putting N electrons per unit cell. 

Now, if originally the system is not at equilibrium and/or if some external perturbation is 
applied, the state of the system evolves in time. Still in the framework of linear empirical models, 
the dynamics is characterized by the unitary propagator Uy{t,s) associated with the effective 
time-dependent Hamiltonian H{t) = -ffpor + v{t): 

7(t) = !7„(t,0)7(0)!7„(t,0)*. 
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Recall that a two-parameter faraily of unitary operators U{t,s) {s,t G R) on L^(]R^) is a uni- 
tary propagator provided (see [17l Section X.12]) (i) V(r, s,t) g R^, U{t,s)U{s,r) = Uit.r); (ii) 
U{t,t) ~ 1 (the identity operator); and (iii) U{t,s) is jointly strongly continuous in t and s. 

Similar considerations hold for mean-field models, although the situation is more complicated 
since both the periodic potential Vpct and the perturbation v depend self-consistently on the 
state 7. The time-evolution corresponding to the Hartree model is known as the time-dependent 
self-consistent field equation, and is in fact equivalent under some assumptions to the so-called 
random phase approximation; see the discussion in |12j . 

We focus in this work on the evolution of the electronic state in insulating (or semiconducting) 
crystals with local defects. The precise functional setting allowing to describe local defects in 
insulating crystals is recalled in Section 12.21 The equation governing the time evolution of the 
defect can be motivated by a formal thermodynamic limit based on the evolution equation ^ 
for finite systems, writing ^(t) = ^'^^^ + Q{t) in the thermodynamic limit, and using the formal 
relation H'^{t) ~ H^^j. + v{t) + Uc(pQ(t))- This leads to the following nonlinear dynamics for a 
given external time-dependent potential v{t): 

i^= [iJp"„+t;e(pQ(0)+^'W,7^er + QW]- (3) 

This paper is organized as follows. After recalling the structure of the time- independent Hartree 
model for perfect crystals and for crystals with local defects in Section [21 we study in Section [3] 
the effective dynamics 

i^=[<, + Mt),7p^ + Q(t)], (4) 

where w(t) is a given effective potential. In particular, we prove that, if the initial condition Q(0) 
belongs to the functional space Q introduced in [S] to describe the electronic structure of local 
defects (see Section [2.2p . and under some reasonable assumptions on the external perturbation w, 
the dynamics is well-posed in Q for all times. We also investigate the linear response corresponding 
to the effective dynamics and show how the results obtained in [5] for the static case can be 
recovered by an adiabatic limit. 

In a second step (Section |4]) , we study the mathematical properties of the nonlinear dynam- 
ics In Scction [4.1[ we prove the global-in-time existence and uniqueness for ^ in the space Q, 
for initial data in Q (corresponding to local defects) . We also provide in Section 14.21 a mathemat- 
ical derivation of the Adler- Wiser formula [1] I23j relating the macroscopic frequency-dependent 
relative permittivity tensor to the microscopic structure of the crystal at the atomic level. This 
derivation is based on a linearized version of the nonlinear dynamics ([3]). Note that a formal 
derivation of the expression of the macroscopic frequency-dependent relative permittivity tensor 
for a general TDDFT dynamics is presented in [TT]. 

The proofs of the results presented in Sections [3] and H] are gathered in Section \5\ 

2 The time-independent Hartree model for crystals 

In this section, we briefly recall the main properties of the time-independent Hartree model for 
perfect crystals and crystals with a localized defect (see [SJE] for a detailed analysis). We consider 
the bulk limit where the nuclear charge of the perfect crystal is described by a 7^-periodic distri- 
bution PpCT, TZ denoting a periodic lattice of M.^. In the sequel, we assume that p^^^ is a locally 
bounded measure. 

2.1 Perfect crystals 

The density matrix 7p|,j. of a perfect crystal obtained in the bulk limit is unique [S]. It is the 
unique solution to the self-consistent equation 

7por ~ 1(-cx3,£f1 (^per)i ^por = ^ ^cr, 
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where Vpcr is a 7?.-periodic function satisfying 

-AFpcr = 47r - p^l'^) , with pI^^.{x) = -f^^^ix, x), 

and where ep € K is the Fermi level. The potential Vpoi is defined up to an additive eonstant; if 
Vper is replaced with Vpa + C, ep has to be replaced with ep + C, in such a way that j^^j. remains 
unchanged. The function Vpei- being in Lp^j.iM.'^), it defines a A-bounded operator on L^(K^) 
with relative bound zero (see [18l Thm XIII. 96]) and therefore H^^j. is self-adjoint on L'^{M.^) with 
domain H'^{M.^). In addition, the spectrum of Hp^^ is purely absolutely continuous, composed of 
bands as stated in [221 Thms 1-2] and [HI Thm XIII. 100]. 

More precisely, denoting by TZ* the reciprocal lattice, by F the unit cell, and by F* the first 
Brillouin zone, it holds 

= U {en,,} , 
n>l, qer* 

where for all q CzT* , {Sn,q)n>i is the non-decreasing sequence formed by the eigenvalues (counted 
with their multiplicities) of the operator 



acting on 



-2 

endowed with the inner product 



Ll^,{r) {u G Ll^{M.') \ u 7e-periodic} , 



= I uv. 



We denote by {un,q)n>i an orthonormal basis of Lp^j.(r) consisting of eigenfunctions of (^Hp^^.j 
The spectral decomposition of {H^^^), thus reads 



^pcr)„ — ' 'y^,^n,q\Un,q){Un,q\- (5) 



' 9 

11=1 



Recall that, according to the Bloch-Floquet theory [TB], any function / e L^(IR'^) can be decom- 
posed as 



where /p» is a notation for |F*|^^ Jp, and where the functions fq are defined by 

Ben ' ' KeTZ* 

Throughout this paper, we use the unitary spatial Fourier transform 

(J-,/)(fc) = 7(fc) = (27r)-3/2 f /(^) e-'^- dx. (7) 



It can be shown that, for almost all q e R^, fq e Lp^j.{r). Moreover, fq+K{x) = fq{x)e '^'^ for 
all K e TZ* and almost all q G R^. Lastly, 



If the crystal possesses N electrons per unit cell, the Fermi level ep is chosen to ensure the 
overall neutrality of the unit cell: 

iV = El{'ZeF* I £„,,<ep}|. (8) 

n>l 
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In the remainder of the paper, we assume that the system is an insulator (or a semi-conductor) 
in the sense that the N^^ band is strictly below the (A^ + 1)*^' band: 

E+ maxEAf., < mineAr+i , S^, 

In this case, one can choose for ep any number in the range (S^, For simplicity we set in 

the following 

2 

and denote by 

9 = S^+i - S+ > (9) 

the band gap. 

2.2 Crystals with local defects 

Before turning to the model for the crystal with a local defect which was introduced in [5], 
let us recall that a bounded linear operator Q on L^(R'^) is said to be trace-class [HI [21] if 
J2i {4>i7VQ*Q4>i)i^2 < for some orthonormal basis of L^(K^). Then Tt{Q) = ^ ■ {(j>i, Q(j>i) 
is well-defined and does not depend on the chosen basis. If Q is not trace-class, it may happen 
that the series ((/>i,Q0i)^2 converges for one specific basis but not for another one. This is the 
case for the operators Q^^ep introduced in ([T5|) (see the results of [5]). 

A compact self-adjoint operator Q = J2i ^i\4'i){4>i\ G S{L^{R^)), with {(pi , ^2 = Sij, is 
trace-class when its eigenvalues are summable: |Ai| < oo. Then the density 

+ 00 

1=1 

is a function of i^(]R'^) independent of the chosen orthonormal basis (0i) and 



+ 00 „ 

TV(0) = ^A,= / 

,;=l 



PQ- 



A Hilbcrt-Schmidt operator Q is a bounded operator such that Q*Q is trace-class. 
We also need to introduce the Coulomb space 



C 



where S^' denotes the space of tempered distributions, the dual of the Schwartz space S^. Endowed 
with the scalar product defined by 



mm,. 



C is a Hilbert space. Recall that L^/^(R'^) ^ C and that, for /i and /2 in 

Difuh)=f I ^^d.dy. (10) 

jr3 f - y\ 

Considering L^(M'^) as a pivot space, the dual space of C is 

C -.^ {V e L^(R3) I WV e (i2(R^))^} , 
endowed with the inner product 

{Vl,V2)c' ■■= ^ f WVi-\/V2^^ f \k\''Viik)V2{k)dk. 

47r Jr3 47r J^s 
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We now describe the results of [Sj dealing with crystals with local defects. The appropriate 
functional space to describe local defects is the convex set 

/C = {Qg Q| -7p°e..<Q<l-7pV}, (11) 

with 

Q = {Q e ©2 I Q* = Q, Q^' e ©i, Q++ e ©i, (12) 
|V|Q e ©2, IVIQ— |V| e ©i, |V|Q++|V| e ©i}, 

where ©i and ©2 denote respectively the spaces of trace-class and Hilbert-Schmidt operators on 
i2(M3) and 

Q- ■■= 7^crQ7°cr, 7°erQ(l " 7^er), 

Q^- ■■= (1 - 7pV)07p°or, Q + + (1 - 7°er)Q(l " 7p°or)- 

Endowed with the norm defined by 

ligils- 11(1- A)i/20||e, + 11(1- A)i/2Q-(l-A)i/2j|e^ + ||(l_A)i/20++(l-A)i/2||g^, (13) 

Q is a Banach space. Although a generic operator Q G Q is not trace-class, it is shown in [5] that it 
can be associated a generalized trace Tro(Q) = Tr((5++)+Tr(Q ) and a density pq € L^(R^)nC 
In addition, the mapping Q 3 Q ^ pq ^ L?iM?) n C is continuous (see Proposition 1]) and 
there exists Cp > such that 

\\pQ\\u^nc<Cp\\Q\\Q, (14) 
for any Q £ Q. Note that if Q e /C n ©1, then of course Tro(Q) = Tr(g), pq G L'^{R^) and 

It is proved in [5] by means of bulk limit arguments that, for insulating and semiconducting 
materials, the ground state density matrix of a crystal containing a local defect, with nuclear 
charge density Pp"J: -I- z^, reads 

7 = 7pcr + Q.,eF- (15) 

The operator Q^^sf obtained by minimizing over /C the energy functional 

E,,,,iQ) = Tro (i/p"„.Q) - / pq{i^*\- T^) + ^^(pq,Pq), (16) 
where Tro (-^pcrQ) a notation for 

Tro (i/p"crQ) = Tr - 6^]'^' (Q++ - Q— ) \H^,, - epr^") + eFTro(Q). (17) 

The energy functional E^,^^ is well-defined on K, for all ly such that (j^* | • |~^) G L'^{M.^) + C . 
The first term of E^^^-^ makes sense as it holds 

ci(l-A)<|i?°,,-eF|<c2(l-A) (18) 

for some constants < ci < C2 < oo (see [5J Lemma 1]). The last two terms of E^^^-^ are also well 
defined since pq G L^(R^)nC for ah Q e /C. 

3 Response to a time-dependent effective potential 

In this section, we study the evolution of the electronic state of the system when the mean-field 
Hamiltonian H^^^ of the perfect crystal is perturbed by a time-dependent effective potential w(i, x), 
so that the system is described by the time-dependent Hamiltonian 

H,{t) = i/p^. + v{t, •) = -^A + 1/pe, + v{t, •). 
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Under the additional assumption 

V e C\R,L°°{M.^)), (19) 

we can apply Theorem X.71 in [17] and obtain the existence of a unitary propagator (Uy {t, io))(io,t)GRxR 
on L2(R^) such that for each ip € H^{R^), and each to e R, t ^ 0to(O Uy{t,to)ijj is in 
C^(R,L^{R^)) n C°{R,H^{R^)), and satisfies 

■ d(t)to (t) , 



dt 



-it) = H4t)cj)toit), (btAto) = ^- 



Besides, denoting by Uo{t) = e~'*'^p<=' the unitary propagator associated with the time-independent 
Hamiltonian H^^j., {Uv{t, to))(to.t)GRxR is the unique unitary propagator satisfying the Dyson equa- 
tion ^ 

V(to,i) e K X R, Uy{t,ta) ^Uo{t-ta)~i Uoit - s)v{s)Uy{s,to) ds. (20) 

J to 

Under the weaker assumption that 

veLUR,L^iW')), (21) 

it can be proved (see Lemma 1151 in Section IS.ip that there exists a unique unitary propagator 
solution to (|20p . By extension, we will call {Uy{t,to))(to.t)eRy<R the unitary operator associated 
with the time-dependent Hamiltonian Hy(t). 

Denoting by 7° the density matrix at time t — 0, wc consider the dynamics of the electronic 
state defined by the evolution equation 

7(t) = f/.(i,0)7°C/„(t,0)*. (22) 

Note that the conditions 7*^ € S{L'^{M.^)) and < 7° < 1 are automatically propagated forward 
in time by (|^ . In addition, if (1 — A)7°(l — A) is a bounded operator, and if v satisfies (fTO]) . 
then (1 — A)7(t)(l — A) is a bounded operator for each t G R, and 7(t) is the unique solution in 
CnK,5(i2(R3))) to the differential equation 

ij(t) = [i/.(t),7(t)], 7(0) =7°. 

Considering v{t) as a perturbation of the time-independent Hamiltonian -ffp^rj it is natural, as 
in the time-independent setting described in Section [2.21 (see in particular the definition (|15p ). to 
introduce 

Q{t) = 7(0 - 70„. 

Using (|20p . (P^ . and the fact that j^^^. is a steady state of the system in the absence of perturbation 
{Uo{t)jp^^.Uo{t)* — 7p(,j.)i ^ simple calculation shows that Q{t) satisfies the integral equation 

ytem+, Qit) = Uoit)Q°Uoitr -i f Uo{t~s)[vis)ntr + Q{s)]Uo{t-srds, (23) 



where (5° = 7° — 7poi.. It is easy to see that under the assumption (pi]) on the effective potential v, 
the above integral equation has a unique solution in C°(R+, iS(i^(M'^))). 



3.1 Well-posedness of the effective dynamics in Q 

We now focus on the interesting and important case when v{t) is the effective potential generated 
by a local defect, that is when 

v{t)=v,{p{t)):=p{t)*\-\-\ (24) 

with p G Ll^^(R, L'^ (R'^) n C). The mapping Vc is an invertible bounded linear operator from C 
to C' , and, according to Lemma [T6l below, it also defines a bounded operator from L^(R^) n C to 
i°°(M3). Hence, ff p e Ll^^{R, L"^ (R^) nC), the potential v defined by ^ satisfies The 
following proposition shows that, in this case, can be considered not only as an integral 
equation on S{L'^{R^^)), but also as an integral equation on the functional space Q. 
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Proposition 1. Consider Q'^ G Q, p €1 Ll^^{M.-^- , L'^ (M.^) n C) and v the effective potential given 
by \24-^ . Then, the integral equation I123\) has a unique solution in C°(K+, Q), and for all t G R+, 
Tro(Q(i)) ^ Tro(0"). In addition, if Q° G IC, then Q{t) e JC for all t S R+. 

The proof of Proposition [1] is based on the foUowing three lemmas. 

Lemma 2. Let Q e Q. Then, for all t e M, Uo{t)QUo{t)* e Q, Tro{Uo{t)QUo{t)*) = Tro(Q), 
and there exists a real constant l3 > 1 ( independent of Q and t ) such that 

^wqwq < \\uomuoitr\\Q < mwa- (25) 

Lemma 3. Let g e L'^{R^)nC and Q e Q. Then, i[vcig),Q] G Q, Tro(i[wc(£'), Q]) = 0, and there 
exists a constant Ccom.Q G IR_|_ (independent of g and Q) such that 

\\i[Vc.{g),Q]\\Q < Ccon.,Q\\g\\L^-nc\\Q\\Q- (26) 

Lemma 4. Let v S C . Then, i[v,^'^^^] G Q, Tro(i[w, 7p(,j.]) = 0, and there exists a constant 
C'com G 1^+ ( independent of v) such that 

||i[".7por]||g < C'com||w||c'. 

The rcsuhs contained in Lemma|3]are estabhshed in the proof of Lemma 5], while the proofs 
of Lemmas [5] and |3] can be read in Section 15.31 

Proof of Proposition]^ As v :~ [p) e L1(M+,L°°(]R3))^ we infer from Lemmas and H that 
the afhne mapping 

Q^-i f ;7o(--.s)[i'c(p(s)),7pcr + Q(s)];7o(--s)*ds 

JO 

is continuous from C°(M+, Q) into itself. The existence and uniqueness of the solution to ([^5]) in 
(7*^(18+, Q) can then be proved by standard techniques (see for instance [H]). The preservation of 
Tro((5(t)) also straightforwardly follows from Lemmas [2j |3] and E) Finally, the fact that — 7pcr — 
Qit) <l- 7°or whenever -7°^^. < g° < 1 - 7°„. can be read off from □ 

3.2 Dyson expansion 

The Dyson expansion consists in writing (formally for the moment) the solution Q{t) of as 
the series expansion 

+00 

Q{t) = Uo{t)Q''u„{ty + QnAt), (27) 

n=l 

where the operators Qn.v(t) are obtained by inserting p7p into (j23p and equating the terms 
involving n occurrences of the potential v. In particular, the linear response is given by 

Oi,.(t)--i/ Uo{t~s)[vis),-f^^^ + Uo{s)Q"Uo{sy]Uoit^s)*ds, (28) 
Jo 

and the following recursion relation holds true: 

Vn > 2, QnAt) = -i / Uoit - «) Hs), Qn-iA^)] Uo{t - s)* ds. (29) 
"'0 

The main result of this section is the following proposition, whose proof can be read in Section [531 
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Proposition 5. Let p £ Ll^^{R+ , {M."^) n C) and v{t) := Vc{p{t)). For each n> 1, the function 
Qn,v defined by i28\) for n ~ 1 and by H29\) for n > 2 is in C''(R+, Q), and, for any n > 1, 
Tro(Qn,u(i)) = for all t G M+. Moreover, there exists a constant C > such that 

Vn>l, ytem+, i|Q„,.ft)||Q </? ^ + '[^°"^ (^C l\\pis)\\L2ncd.^ , (30) 

and the right-hand side of {21^ converges in Q, uniformly on any compact subset o/M+, to the 
unique solution to (23^- in C°{R+, Q). 



It is possible, and convenient for some calculations, to reformulate the dynamics (|23p in the 
so-called interaction picture (the reference time being fixed to to = 0), introducing the operators 

f/i„t(t) = C/o(t)*C/,(i,0) and «;„*(*) = C/o(t)*f(t)C/o(t). (31) 

The Dyson expansion of the evolution operator C/int(0 then reads, in terms of the potential in the 
interaction picture, as 

Uint{t) = 1 - i / Vint{s)Uint{s) ds 

Jo 

= l + V(-i)"/ / .../ z'i„t(ii)vint(i2)...Wint(i„)di„...dti. (32) 
Jo Jo Jo 

Note that, in the last integral, the times are increasing from the right to the left (t„ < t„_i < 
• • • < ^i), and the operators ('yint(^j))i<j<n do not commute. We can also rewrite the recursion 
in a form reminiscent of the Baker-Campbell-Hausdorff formula: for any n > 1, it holds 



Qn,vit) = (-i)"f^0(t) { I [«i„t(il), bint(t2), ■ • ■ , ^int (in) , 7per + Q"] ■ • ■ ]] dti . . . dt^ Uoit)* . 

\Jo<t„<-<ti<t J 

3.3 Linear response and definition of the polarization 

The aim of this section is to motivate, using rigorous mathematical arguments, the formula (j44p 
for the polarization matrix usually encountered in the physics literature, known as the Adler- 
Wiser formula PQI13] (up to a factor 2 accounting for the spin, see (2.8) in [T]). These expressions 
are established for a modified linear response involving some damping. Proposition [S] gives a 
mathematical meaning to the polarization formula when the damping vanishes. We therefore 
focus on the linear response term, which is the operator Qi,t,(t) given by ([25)) : 



gi ,„(t) = -i / U^{t - s) {v{s),-il^^ + t/o(s)Q"f/o(.s)*] f/o(i - s)* ds. 





We choose = 0. When the external perturbation v[t) is compactly supported in time in some 
interval [— ioj^oli we can view the perturbation process as a dynamics starting in the distant past 
from an equilibrium state described by Qify = up to time t = —to, and perturbed only for times 
t > —to- Upon changing the reference time from to —to, the following integral equation is then 
obtained: ^ 

VieR, Qi,„(i)=-i/" Uo{t-s)[v{s),^^^^]Uo{t-s)*ds. (33) 

The interest of this formulation (compared to the original formulation ((28|) ) is that it can be 
interpreted as some time convolution, which can then be rewritten in a simpler manner using 
Fourier transforms in time. Using Lemmas [5] and HI and the density of C^(R,C') in L^(R,C'), it 
is easily seen that v defines a linear mapping from i^(M,C') to C{J(R, Q), where Cj](R, Q) 

denotes the space of the continuous bounded Q-valued functions on R. It is then possible, by 
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density, to consider external perturbations v £ L^{M.,C'), and not only perturbations with compact 
supports in time. Alternatively, for a given perturbation v[t) defined only for positive times, the 
linear response can be written as p3|) upon considering v{t) = v{t) if i > 0, and otherwise. 

Since Qi.v{t) G Q for all t e M, it is possible, in view of [SJ Proposition 1], to associate a 
density PQ^^^t) e L'^{R^)nC to this operator. This defines a bounded linear mapping 

Xo ■■ L\R,C') C^^{M.,L^{m.^)r\C) 

1/2 1/2 1/2 

In fact, it is more convenient to work with the mapping S' = Vc XoVc ■ As Vc is an invertible 
bounded hnear operator from L^(R^) onto C, and from C onto L^(R^), and as L^(R^) n C = 
H'^iM.^), is a continuous linear operator from Li(R, 2.2(^3-)-) C0(R,7Ji( 



f ^ 1/2 / (34) 



In order to state our results, we need to introduce additional Fourier transforms, taking the time 
variable into account. The partial Fourier transform with respect to the time variable, denoted 
by Ttf , lias the following normalization^ 

/ + 00 
f{t,x)e'^'dt. (35) 
-oo 

The space-time Fourier transform J-t^x based on J-t and on the spatial Fourier transform J-^ defined 
in is then 

(J-t,,/)(c.,/c) = (J-tJ-,/)(a;,fc) = (2^)-3/2 /■ /(t,x)e-'('=-^-'^*)didx. (36) 



3.3.1 Damped linear response 

In order to study the properties of the linear response, it is convenient to first focus on the damped 
linear response defined, for > 0, as 

QUt) = -i /* Uoit ~ s) [v{s), 7^er] Uo{t s)*c-"(*-^) ds. (37) 
We denote the associated damped linear response operator 

f ^ 1/2 \ (38) 

As shown below (sec Proposition[5]), the operator S"^ indeed is an approximation of the operator S". 
The interest of the operator is that it has better regularity properties than the plain linear 
response S". 

For a given 77 > 0, we consider a simple closed contour in the complex plane, symmetric 
with respect to the real axis, enclosing a(Hp^j.) n (—00, ep], containing no element of K. ± h] (see 
Figured]), and such that 

dist('4,a n (-00, ep]) > |, dist(^^,R + ir?) > |. (39) 



We can then prove the following result. 

^Note that, as usual in the physics literature, there is no minus sign in the phase factor in the definition of the 
Fourier transform with respect to the time variable. 
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177 




Figure 1: Integration contour used in Proposition [HI 



Proposition 6. The operator is bounded on L'^{R, L'^{R^)) and, for /i, /2 e L'^{R, L^{M.^)), 

(/2,^Vi)i.(i.)= / (-Ft/2M,^"M^t/iM)i2(K3)du;, (40) 
where, for /ii, G L^(R'^), 



(/ii,(f''(a;)/i2)i2(R3) = -Im 



Tr 



1/2/1 N '''per 1/2 /J, N 



dz 



In addition, there exists a constant C £ IR.+ sueh that 

II^''IIb(L2(RX2(R3))) = sup ||(?''(w)||e(i2(K3)) 



(41) 



< 



c 



The proof of this resuh can be read in Section [S75l The operator appearing in the trace on the 
right-hand side of (PT|) is indeed trace-class since 



Tr 



(7p"e..)^ 



per 



Tr 



L2(R3) 



(7°e.)- 



^ - (^^cr + ^ + i?7) 



1/2/1. \ 
I'c (/i2),7po 



7° 
/pc 



z- H° 



7pV,«y'(/^i) 



(42) 



and [vj (^),7por] ^ ®2 when h G i^(M'^) by LemmalU In addition, in view of the conditions p9p 
on the contour "^^j the operators ^_f^jj^'°^j^^j^-^^-^ and are bounded uniformly in z and uj, 

with a bound proportional to 77"^ for both of them. The right-hand side of (PO)) is therefore well 
defined for /i,/2 G ^^(m, ^a^j^S)) 

Since the linear response commutes with time translations, it is not surprising that the op- 
erator S'^ is diagonal in the frequency domain (in the sense of (001) )■ Moreover, the operators 
£"^{00) commute with spatial translations of the lattice. They are hence decomposed by the Bloch 
transform associated with the lattice TZ. The action of (o'''(a;) on the fiber associated with the 
Bloch vector g € F* is denoted by S'^{uj,q). Introducing the Fourier basis (e_R-)_R-g7^. of Lppi.(F), 
where exix) = |F|~^/^e''^ '^, the Bloch matrices of the operator S'^'(u)) are defined as 

and it holds 

VA-e7^^ J-t,, (/?"/) ('^, 9 + i^) - J2 '^K,K'i^^l):Ft,Jiuj,q + K'), fora.a. (w,g) gMxF*. 

(43) 

As stated in the proposition below, the Bloch matrices of the operators (o^ can be written 
in terms of the Bloch decomposition of the mean- field Hamiltonian Ifpgr- The corresponding 
expressions are known in physics under the name of Adler- Wiser formula [2 123) . 
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Proposition 7. For each rj > 0, the Bloch matrices of the damped linear response operator S'^ 
are given by 

( \ lr»(g) \q + K'\ 

'^K.K'K^^l)^ \qj^K\ ^K,K'i^'l)^ (44) 

where the continuous functions ^, : M x E'^ — C defined by 

^KK'\^il)~ y. {^n<N<m^i-m<N<n) 'f : dq 

— , Jr* £n,q+q' — £m,q' — LO — \1J 

n,m—l i-i I y '-i 

(45) 

are uniformly bounded. 



The proof of Proposition [7] can be read in Section 



The above expressions make sense 



since it is proven in Lemma \W\ that the sums over m, n which enter in the definition of S'j^ j^, are 
convergent. This is due to the fact that for all r] > 0, (i) en,q grows as n^/'^ when n goes to infinity, 
uniformly in q G F* (see ((78)) ): and (ii) for a given K S TZ* , there exists a constant C > such 
that, for all 1 < n < TO > AT + 1, and q, q' e T* (see pi)) ). 



iK-. 



^n,q+q' /i^2 



<Cil + \K\^) 



For later purposes, it is useful to notice that for all / G S^{M. x M.^), 

:Ft.ASV)= ^K{r^K'{J'u.f)<^lK'), (46) 

where Txfi^, q) — f{^^, q — K) is the momentum space translation of vector K. As will be shown 
below (see Proposition [5]), this representation is well suited to the limiting procedure 77 — ?> 0. Note 
that (f^ j^, belongs to L°°(Mx M^) and hence defines a tempered distribution on Rx M^. Therefore, 
T-K' {p^t,xf) K' ^ tempered distribution when / e .5^(R x R^). The fact that the series on 
the right hand side of converges to Ft^x (S'^f) in the sense of the tempered distributions is 
proved in Lemma [23l 

The expression ([15)) is a result of the following computations. Since the (^'jj^^/'s are C- valued 
functions on M x M"^ with supports in M x F* and uniform bounds in L°°{R x M.^) (for fixed 77 > 0), 
we obtain that, for aU / e J5^(R x R^) and all ip e ^ (R x M.^) , 



{:Ft,x {^V) .H>)sj"sy= / (^V) (^, fc) ¥'(w, fc) rfw dk 

^ f J2 f ^t,xi'^y)iuJ,q + K)ip{LU,q + K)dqduj 

= f f SlK,{u:,q)Ft,xf{^.q + K')^{u,q + K)dqdu: 

^ Y j^,T_K'{J't,xf)T-KV 



.9" ,.9 
K.K'en* 



Y TK{T-K'{J't,xf)'S'l^K')^'p) 
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3.3.2 Bloch matrices of the linear response 

In order to characterize the Bloch- frequency decomposition of the operator S defined by (|34p , we 
investigate in this section the limit of the damped linear response when ?/ J, 0, by passing to the 
limit in 

Proposition 8. The operators S"^ converge to S in the following sense: for any f G _L^(]R, i^(R'^)), 
VTeM, lim = in L°°{{-oo,TIH\R^)). 

In addition, for each {K,K') G TZ* x TZ* , the Bloch matrix S'^ j^, converges in 5^'{M. x M"^), when 
rj — >■ 0, to a limiting distribution denoted by Sk.k' ■ Finally, for each f £ ^(M x R"^), the following 
equality holds in .y"{R xR^): 

K' {^t,xf) Sk,k']- (47) 

The proof can be read in Section [5771 This result shows that the matrix {(^k,K')k.K' can be 
interpreted as the Bloch matrix of the operator S". An expression of {<^k.k' t'^) jyi y is provided 
in the proof of Lemma [H] 

A careful inspection of the proof shows that (|47|) can be given a meaning for functions / which 
are not in ^(R x R"^), but are nevertheless regular in space and decaying in time, see Remark [Ml 



Remark 9. The tempered distribution Sk.K', defined in Proposition\^ as the limit of S"^ when 
rj goes to zero, can be written more explicitly when the pulsation ui is not too large, namely when 
its absolute value is smaller than the band gap g defined by Indeed, when \u)\ < g, it holds 

^ ^\ ^ 9 ^ \^\ > ^ for all q,q' and all n,m satisfying l<n<N<mor 



'-■n,q+q' 



1 < m < N < n, so that for all K, K' G TZ* and almost all (cj, g) G M x F; 

Ir-(g) \q + K'\ ^0 

Sk,K' (w, q) = 

where the bounded continuous functions ^/ : M x M"^ — > C are defined by 

njn—l 



(48) 



^n.q-{-q' ^rti.q' 



dq'. 



Let us also notice that 



\q + K' 



-S'K,K'i(^,q) 



Ir-(g) 
27ri 

lr»(g) 
27ri 



Tr 



-iKx 



(7pcr)g+9' 



AK'x 



(49) 



(7° )^ 

dq'dz 



Tr, 



-IT 



(7p°e.),' 



(^pV) 



9+9' 



-(^p°e 



- o?g' dz 
(50) 

where 'lo^ is a contour enclosing {a{H^^^) n (— oo, ep]) ztoj and containing no element of a (iJp^r) 
[eF,+oo) (see (|83p below and Figure\^. 



3.3.3 Adiabatic limit 

The linear response of the electronic density for time- independent perturbations was studied in [5] . 
The aim of this section is to recover the static polarization operator from the time-dependent one 
in some adiabatic limit. 
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Figure 2: Integration contour used in formula ([5n|. 

The static linear response operator is defined as := v^^'^y^Q^'^'^'^v\^'^ ^ where is the 

time- independent polarizability operator introduced in [21 Proposition 1] : 

with (for F e C) 



1/2 / A (51) 



Qstatic ^ 1 i i^p^)- V(Z - dz, 

where "^^o is a simple closed contour in the complex plane enclosing a(il^^^ H (— oo,£f] and con- 
taining no clement of a (-ffpor) ^ I^f, +oo). Wc deduce from and the results in [5] that 

'i{K, K') e 7^* X 7^^ S^K'^K'i^) = ^K.K'(0, g) for a.a. g € T*. 

The timc-indcpcndcnt polarizability operator is therefore the zero-frequency limit of the dynamical 
response. 

This consideration leads us to study the adiabatic limit of the linear response. To this end, we 
consider the following time evolution for some parameter a > small enough: 

QIM = -^I^ t^o (^-s) [«(«s),70J C/o (^^-5^ ds. 

In the above dynamics, the evolution of the time-dependent potential v is slowed down, and the 
effect of the perturbation is considered on longer times (for t > 0) in order to obtain a non- 
trivial result (note that Qf y{t) = Qi.v(a-)it/ct), where Qi^v is defined in l\'S'3^ ). Equivalently, this 
procedure may be seen as accelerating the free evolution generated by H^^^ and appropriately 
rescaling the result. Indeed, a change of variables shows that 

For any a > 0, we introduce the rescaled linear response operator 

: L1(R,L2(k3)) ^ C{;(M,i7i 



'(/) 



Proposition 10. For any function f e .y(R x R^), 

lim(?"/ = (?"/ iny'i 

where for all t G R, {S^ f)(t) = S^^'"^'" {f {t)) . 

This result is proved in Section 15.81 It means that, in the adiabatic limit, the linear response at 
time t only depends on v{t). There is no memory effect. Moreover, the linear response at time t is 
given by the time- independent (or static) polarization operator S'^^'^^^'^ studied in [SI Proposition 4] . 
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4 Nonlinear Hartree dynamics 



We now focus on the nonUnear Hartree dynamics defined by 

yteW+, Q{t) = UoWUo{ty-i [ ?7o(t-.s)k(PQ(.)-Ks)),7pcr + Q(.'*)lf^o(t-s)*ds, (53) 

Jo L J 

for an initial condition Q° G /C, and for a nuclear charge distribution of defects i'{t) e L'^{M.'^) n C 
for all t. Recall that the solutions of ([5^ are the mild solutions of the von Neumann equation ([3]) 
with v{t) = -Vciiyit)), since [i?°er>7por] = 0. 

4.1 Well-posedness of the dynamics 

The main result of this section is the following. 

Theorem 11. LetuG Ll^^{R+,L'^{R^))r]W^^^{R+X)- Then, for anyQ° e K., the time- dependent 
Hartree equation (j53p has a unique solution in C'^(M+, Q). Moreover, for all t > 0, Q(t) G K. and 
Tro(Q(t)) = Tio{Q'^). Finally, j/7por + "■''^ orthogonal projector, then 7°^^ + Q{t) is also an 

orthogonal projector for all t > 0. 

The proof of local existence and uniqueness (see Section I5.9.1|) is classical and is based upon 
a Banach fixed point argument in a well-chosen ball of Q. Once local-in-time existence and 
uniqueness is ensured, it is possible to extend the well-posedness of the dynamics to all times by 
proving that the Q-norm of Q{t) does not blow up in finite time (see Section [5.9.21) . This can be 
performed by controlling the growth of || (5(t) || q by means of the energy functional £ : M+ x Q K 
defined by 

£{t,Q) = i?.(t),ep(Q) = Tro(i?°,,Q) - D{pQ,,,it)) + ^D{pq,pq), (54) 

where Tro{Hp^j.Q) are defined respectively in ((T5)) and pT)) . 

Under appropriate regularity assumptions on Q'^ and v, the unique solution of (j53|) is a classical 
solution of ^ with v{t) — — Uc(i^(i)). Let us detail this point. The evolution problem ([55]) can be 
formally written as 

&l=AQ{t) + F{t,Qit)), Q(0)=QO, (55) 

where 

F(<,Q)=-i[t;e(pQ-i^(t)),7pcr + Q], 

and where A is the generator of the strongly continuous group (G(i))tgR on B{Q), the space of 
the bounded linear operators on Q, defined as 

G{t)Q^U'\t)QU'>ity. 

In fact (see [TOl Section XVII.B.5.1]), 

D{A) = {q e Q QD{H°J C D{H^J, -i(ffpcrQ - Q^pcr) is a linear operator on ^2(^3)^ 

with domain D(Hp^^.), which can be extended to a bounded operator on L'^{M^ 
and this extension, denoted by Qa, belongs to q|, 

and AQ — Qa- In particular, the operators Q G Q such that (1 — A)Q(1 — A) e Q are in D{A) 
(recall that the set of those operators is dense in Q, see [SJ Lemma 2]), and AQ — —i[Hp^^., Q] in 
this case. We can then state the following. 
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Proposition 12. Let u G C^(M+, ^^(M^) nC). Then, for any Q° e D{A), the unique solution to 
the time- dependent Hartree equation ()53p in C*'(K+, Q) is in C'^(M+, _D(A)) H C"'"(R+, Q), and is 
a classical solution to I155\) . 



Proof. The result follows directly from Theorem 6.1.5 of [TB], since F : M+ x Q — > Q is C^. 
Indeed, this mapping is differentiable and its derivative 

dF{t,Q) ■ {S,R) = isK(i''(i)),7per + Q]- ibc(pQ " ^^W),^] - i^c (pi?) , 7pcr + Q] 

defines a bounded linear operator on ]R+ x Q. Moreover, in view of Lemma [31 the mapping 
(t, Q) dF{t, Q) is continuous whenever v e C^(R+, L'^{R^) n C). □ 



4.2 Macroscopic dielectric permittivity 

We start with formal computations, which, to be justified, would require estimates on the long 
time behavior of Q{t). Unfortunately, we do not have such estimates, see the discussion after 
Proposition [T31 For the same reasons as the ones presented before ([M]). we choose = in (|53p 
and change the reference time from to to, letting then to go to — oo, formally obtaining 

Q{t) = -\f Uo{t-s)\v,{pQ(^s)~v{s)),ll.r + Q{s)\Mt-srds. (56) 

The above integral equation can be rewritten as 

Q{i) = Qi,v^(pQ-v){t) + Q2,v^(pQ-u){t), (57) 

where the linear operator (5i,u is defined in (j28p . and where the remainder Q2,v^{pQ-v){t) collects 
the higher order terms. Equation (|57p can be reformulated in terms of electronic densities as 



PQ{t) = [Ciiy-pQ)] it)+r2it), 
where C = — Xo'^c and r2(t) — or equivalently as 

[{l+C){u-pQ)]{t)^:.{t)-r2{t). (58) 
This motivates the following result (proved in Section [5.10|) . 



Proposition 13 (Properties of the operator £). For any < ft < g (the band gap of the host 
crystal), the operator C is a non-negative bounded self-adjoint operator on the Hilbert space 

^2 = i'(K,C)|supp(J■t,,^)) c [-0,17] x K^j, 

endowed with the scalar product 

{Q2, Qi)l^(c) = / D(^g2{t, •), Qi{t, •)) dt = A-K 



Ft,xQ2{t^, k)Ft,xgiiuj, k) 



duj dk. 



Hence, 1 + £, considered as an operator on Ji(ti, is invertible. 



This result cannot be used as such to study ([55)) since, even when belongs to J^n {0 < V, < g), 
the nonlinear response r2 generally involves frequencies with absolute values larger than 12. This 
can be seen from the relation (|29p . For instance, Q2,v is a convolution between the time evolution 
Uq of the perfect crystal, and products such as vQi^y. Since the time Fourier transform of each 
of the element of the latter product has support in (—17,17), the time Fourier transform of their 
product has support in (—217,217). 
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In order to rigorously obtain the macroscopic dielectric permittivity from (j58p . some spatial 
rescaling should be performed. In the time-independent case dealt with in [6], the equivalent of 
the nonlinear term r2 turns out to become negligible under this spatial rescaling. In order to 
prove that the same phenomenon occurs in the time-dependent case, we would need estimates 
on the time growth of the nonlinear term r2{t). Controlling this term is probably difficult since 
very few is known about the long time limit of dynamics such as (j56p . Typical tools to this 
end are Strichartz-like estimates, which allow to establish appropriate decays in time and prove 
scattering results (see for instance [19l Section XI. 13]). Such inequalities are easy to prove for 
the operator — A on L^(M'^). To our knowledge, the only known dispersion inequality for periodic 
Schrodinger operators is restricted to the one-dimensional setting, see the recent work [3]. The 
proof is based on the stationary phase method, but several fine estimates rely explicitly on the 
fact that the system is one-dimensional. It is unclear whether such results can be extended to 
three-dimensional systems. 

We will therefore limit ourselves to pass to the macroscopic limit on the following linear prob- 
lem, obtained by neglecting r2 in ((58| : < < g and u <E 'Hn being given, seek e Ho such 
that 

VteR, [{l + C){v- p^)]{t) = v{t). (59) 

In order to study the response of the system at the macroscopic scale, we consider the regime 
where the perturbation is weak but spread out over a large region, using the same spatial rescaling 
as in [5]. For > 0, introduce the rescaled charge of the external perturbation 

Uri{t,x) = rfv{t,rix). (60) 

Note that i^rjit, x) dx = J^3 v{t, x) dx for all 77 > and all i S M. We also define the rescaled 
potential generated by the total charge of the defect — p^^ as 

w^;'(t,x) = r^v,{vr,^p,^) {t.r^x) . (ei) 

The scaling of the potential is such that, in the absence of dielectric response [C = 0), the potential 
effectively seen by the crystal is WU = Wc(j^)- We are then able to prove the following result. 

Proposition 14. There exists a smooth mapping {—g,g) 9 w i— >■ eM('^); with values in the space 
of symmetric 3x3 matrices, satisfying eM(w) > 1 for all lu € {—g,g), such that, for all v g 'Ho 
(D < 17 < gj, the rescaled potential defined by ([SHI - lpT]) converges weakly in T-Lq when r/ goes 
to 0, to the unique solution in Jifi to the equation 

- div(eM(w)V [J^tW,] (uj, •)) = 4^ [-^tH (oJ, (62) 

where div and V are the usual divergence and gradient operators with respect to the space vari- 
able X, and where Tt is the time Fourier transform defined in (j35p . 

This result is proved in Section fS.llI In particular, the precise expression of eM(w) in terms of 
the Bloch decomposition of the mean-field Hamiltonian H^^j. is given in (|110p . Note that in the 
macroscopic equation the pulsation w enters as a parameter: there is no coupling between 
different values of uj. In the space-time domain, this means that the charge ;/(t, x) and the potential 
W,y(t,x) are related by a space-time convolution. 

5 Proof of the results 
5.1 Existence of propagators 

Lemma 15. Consider a self-adjoint operator H^^j., and the associated propagator Uo(t) = e^'*^p=' . 
For a given potential v G LJq^(R, L°°(R'^)), there exists a unique unitary propagator {Uy(t, s))(^s,t)£«2 
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satisfying the integral equation 

V(t,to)eMxK, Uy{t,to) = UQ{t -to) -i [ Uo{t~ s)v{s)Uy{s,to)ds. (63) 

J to 



Proof. The proof uses ideas from [ITj Section X.12], [TBI Section 4.1] and [24]. Denote Wint(i) = 
Uo{t)* v{t)UQ{t) the perturbing potential in the interaction representation. Note that vint belongs 
to Ll^^(M, 6(L2(M3))), where B{L'^(R.^)) is the Banach space of the bounded operators on L'^(R^). 
We first define the propagator f7„,int(i, ^o) associated with the family of bounded operators Wint- 
In view of |T71 Section X.12], we can associate to the family of bounded operators (wint(O)teR a 
unitary propagator (?7„,int(i, io))(to,t)eR2- In addition, 

y{to,t) G M X R, Uv^int{t,to) = 1 — i / Vint{s)Uv^int{s,to ) ds. (64) 

J to 

It easily follows that (t/t,(i, io))(to,t)eR=: where Uv{t,s) ~ Un{t)Ui,^int{t,tn)Uo{to)* , forms a unitary 
propagator, and verifies ([55)1 in view of (|64p . □ 

5.2 Some properties of the Coulomb potential 

The following result is an extension of Lemma 3 in |13j . 

Lemma 16. When g e L'^{R^) n C, the potential Vg = Vc{g) belongs to L^{R^) n L°°(M^), and 
there exists a constant Cpot G K+ such that 

Vpe[6,+(»], \\Vg\\Lp <Cpot\\g\\L-nc- (65) 

Moreover, for all q G [2,6), VVg € ^L^{R^) n L''(]R^)^ , and there exists a constant Cg^M,q G ^+ 
such that 

||VVg||L5 < Cgrad.qllellL^nC- 



Proof Note first that, since Vg G C, it holds Vg G L6(M3) and Wg G (L2(m3))3 ^j^j^ 

WVghe < Cell g\\c, m\\L-^<Cc\\g\\c, 

for some constant Cq G M+ independent of g. The boundedncss of the potential comes from the 
following estimate: 

V 71- J|fe|<_R l^r V TT j\k\>R Fr 

Now, 
and 

/ i^*^f/ ^^^V'^yf ii^ii-' 
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so that finally 



2R 

VsWl^ < \l j^\\Q\\L^ + \l —\\g\\c 



By interpolation, Vg G L^'(R'^) for any p e [6, +00], and the constant Cpot in ((65)) can be chosen 
independently of p. 

To show that VVg S (L'?(M3))3 fo^. any 2 < q < 6, it is sufficient, by the Hausdorff- Young 
theorem (see for instance [l7]), to verify that k ^ F{\W\Vg){k) = |fc|Vg(fc) = 47r^(fc)/|fc| is in 
L"(M3) for 6/5 < a < 2 since 



IVIFglU, <47r(27r 



)3/2-3/q 



\k\ 



with a ^ = 1 — q ^ . Let i? > 0. First, the Holder inequality (with exponent 2/ a and conjugated 
exponent 2/(2 — a)) ensures that 



\m\Y 



\k\<R. \ \k\ 
Second, 



\k\>R 



dk < 



m)\Y 



k\< 



a/2 



dk 



l-a/2 



dk 



< (47r)^-" — 



|fc|<-R 



3\ l-"/2 



\k\ 



dk < 



1/2 



k\>R 



\g{k)\^ dk 



I-q/2 



|fc|>_R 



|fc|2a/(2 



■dk 



< 47r- 



2-a 
6 — 5q! 



i-q/2 



^(6-5«)/2||^||a 



L2- 



when 2a/ (2 — a) > 3, i.e. a > 6/5. 



□ 



We will need the following result. 

Lemma 17. Let V G L°°{M.^) be such that W e {L'^{M.^) f. Then (1 - A)i/2y(i ^ A)-i/2 ^ 
bounded operator on L^(M^), and there exists a constant C e M+ independent ofV such that 



|VV||l4). 



(1 - A)i/V(l - A)-i/2 < C (lll^lli. 
/n particular, there exists a constant C'a G IR+ suc/i f/iaf 

ygeL\R^)nC, (l-A)i/2„^(g)(i„A)-i/2 < CA||^^|U2^c. 



(66) 



Proof. First, note that, for a given smooth function 7/;, 

9,.(i-(i- A)-i/v) = (a,,y)((i- A)-i/2v') +ya,.(i- A)-i/2^. 

The operator (9j;;y)(l — A)~^/^ is in 64 by the Kato-Seiler-Simon inequality (see [101 HI]): and 
is therefore bounded. The operator Vdxi{l — A)~^/^ is clearly bounded since V € L°°{R^). 
Therefore, the operator dxiV{l — A)~^/^ is bounded for i G {1,2,3}. More precisely, there exists 
a constant C G such that the bounded operator A = V{1 — A)~^/^ verifies 

\\A\\ < C\\Vh^, \\dx^A\\ < c{\\Vh^ + \\VVh^). (67) 
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Then, for a given function -0 G L (M'') 



L2 



L2 



L2 



2 

L2 



□ 



L2 ^ II 

i=l 

i=l 

which, in view of (|S7)) and Lemma [TBI gives the expected results. 
5.3 Some stability results 

Before providing the proofs of Lemmas [21 [31 and [TSl we first show that, for any fixed c < 
min(cr(i/pj,j.)), the norm 

(68) 

is equivalent to the norm ||Q||g defined in p^ . More precisely, 

b''\Q\Q<\\Q\\Q<b\Q\Q, 
with 6 = max(||S||, WB-^'^) > 1 where 

i3 = (F°,,-c)V2(f-A)-V2 
is bounded and invertible. This is a consequence of the following inequalities: 



(69) 



and 



as well as 



(iJ° -c)i/2g||g^^||B(l-A)V2g 



(i-A)^/^g = i?-i(ij,L-c)^/^g 



< \\B\\ 



< \\B- 



to2 



(1-A)1/2q 



62 



{Ht,-cf/^Q^^{Ht,-cY/^ ^ = i?(l-A)i/2Q±±(i_A)i/2B 



tol 



Si 



< \\Bf (1- A)i/2g±±(i_ A)i/2 



and 



(l-A)i/2g±±(i_A)i/2 ^ <||i3-i||2 .-c)i/2g±±(i/"„,-c)i/2 



Si 



It is therefore sufficient to prove the stability results we need in the norm | • | g defined in . The 
interest of this norm is that it simplifies some algebraic computations since any function of -ffper 
commutes with {H^^^ — c)^/^. 

The first stability result, stated in Lemma [2l shows that the space Q is stable under the action 
of the propagator of the corresponding periodic mean-field Hamiltonian. 



20 



Proof of Lemma [H The inequality (j25p is a straightforward consequence of the equivalence of 
norms (|68p and the equality 

\uomuo{tr\Q^\Q\Q. 

Moreover, as j*^^^. and UQ{t) commute, we obtain 

TTo{Uo{t)QUo{tr) = Tro(([/o(t)gC/oW*)"") +Tro((C/o(t)Qt/o(t)*)++) 

= Tr(7pVC/o(t)Qf/oW*7p°cr) + Tr((l - jl,)Uo{t)QUo{tni ~ 7p°cr)) 
= TT{Uo{t)Q-'Uo{ty) + TiiUo{t)Q++Uo{tr) 
-Tr(Q— )+Tr(Q++) =Tro(Q), 

which completes the proof of Lemma [51 □ 

The second stability result, stated in Lemma[3J shows that for all g G L^(R'^)nC, Q i~> i[wc(p), Q] 
defines a bounded linear operator on Q. 

Proof of Lemma\^ First, since Vg := Vc{q) is bounded in view of Lemma [T51 and Q € 62, i[Vg, Q] 
is self-adjoint and Hilbert-Schmidt on with 

||i[^„Q]|le, <2Cpot||g||L2ncl|0||s. 

In addition, 

where (i/O^, - c)1/2q ^ ©3 and 

= [Hi, cY'^Vg{Hl, c)-i/2 = 5(1 - A)VV,(1 - A)-!/^^-! (70) 

(with B defined in ([511)) is bounded by Lemma HT] This shows that {H^, - c)i/2i[ye, Q] € ©2 
with 

(i/0,,-c)i/2i[y„Q] < {C^ot + \\B\\\\B-^\\C^)\\gU2r.c\QU 
S2 

Now, consider for instance (i[Vg,(5])++ = i(l — ^'^^,)\Vg,Q\{l — 7per)- The goal is to prove that 
(ijO„ - cY'^iAVe. Q])++(i?o„ - c)i/2 e 61. This operator can be decomposed as 

- \{Hl, c)V2(i _ 7^„,)Q70,,l^,(l - 7^er)«r - 
+ i(i/°,, - c)V^(l - 7^,,)y,Q++«, - cf/' 

- ml. - c)^/^g++F,(i - ji^mi, - cfi\ 

Let us deal with the first and the third terms on the right-hand side (the second and the fourth 
terms are the adjoints of the first and third terms respectively). It holds 

i«, - cfl\\ 7p°cr)^.7p°crQ(l - 7p°c..)«.. c)^'^ = cf'^^Vg, 7^„.] Q'^ {H% ~ cf'^ 

and 

[Hi, - cf'\l - 7°,,)^,Q++«, - of" = (1 - 7p°c.>«. - c)^/^g++(ff°„ - of" 

with w defined in (|70p . In view of Lemmas |3] and [T71 we infer that the above operators are 
trace-class and that 

||i(J?°3, - cf'\l - 7p°cr)^e7°crQ(l - 7;!er)(i^pcr " cf'^\\e^ < Ccom||Sl| 1 1 | L^nC I Q I Q 

||i(ffp"„ - c)i/2(i „ 70„)V^^g++«r - c)i/2||6, < CaIISII ||^^|U^^clQlQ• 
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Using similar manipulations for the other terms, we finally obtain ()26[). Besides, 

(i[14, Q])- + Q])++ = i {V--Q- - Q-V--) + i {V-+Q+- - Q-+v;-) 

It follows from the cyclicity of the trace that Tro(i[Vg, Q]) — 0. □ 

The last lemma of this section is concerned with the regularization operators 

Rs = (l + 6\H"p,,-eF\r'. 

The properties of these operators we will make use of are collected in the following lemma. As 
in [5], we introduce the space 

6? := {Qe&2\ Q++ e 6i, Q-- e ©i} , 

and denote by Tro{Q) = Tr((5++) + Tr(Q^^) the generalized trace of an operator Q € &i. 
Lemma 18. The regularization operators have the following properties. 

(1) For all Q ^ Q and all (5 > 0, RsQRs G Q, and there exists a constant C independent of Q 
and S such that 

yQeQ, y5>0, WRsQRsWq < C'WQWq. 



In addition. 



\im\\RsQRs-Q\\Q = 0. 

6^.0 



(71) 



(2) For all Q £ Q and S > 0, i[Q,Rs] G Q, and there exists a constant C independent of Q and S 
such that 

yQeQ, V<5>0, \\i[Q,Rs]\\Q<C\\Q\\Q. 

Moreover, 

]im\\i[Q,Rs]\\Q = 0. (72) 

(3) Let g e L^R^) C, Vg v^{g) and Q e Q. Then for all 6 > 0, VgRsH^^^^Q, VgQRsH^^^^, 
VgRsHp^j.QRs and VgRsQRsHp^j. belong to &i, and the following estimates hold, for a con- 
stant C independent of g, Q and S: 



Tro(y,M°,,gj <C||e|U.ncllQllQ, Ti'o [VgQ Rs H^^^^j < CWgh^ncWQWQ, 
Tvo(VgRsH^,,QR5)\ < Cllf^llL^ncllQllQ, \Tro{VgRsQRsH^,,) \ < CllellL^ncllQlls- 



Proof. We prove the bounds in the norm defined in 



Let Q E Q and S > 0. It is clear that 



RsQRs is Hilbert-Schmidt and self-adjoint. In addition, {RsQR, 



RsQ Rs- Using the fact 



that Rs is a bounded self-adjoint operator satisfying < i?^ < 1 and commuting with -ffpon 
obtain 



- cy^'{RsQRs)^^{H^^,, c)V^||e^ < \\(H^^^^. c 



per 



c)^/^||e. 



Likewise, 



UH^,, - cy/'{RsQRs)\\l, = Tr(Rs{Hl, - cY''QRIQ{HI, - cY'^Rs) < ||(i?p^ - cy/'Q\\ 



lis 



Hence, RsQRs & Q and \RsQRs\q < \Q\q- The property ([7T|) is established in the proof of [5l 
Lemma 2]. 
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Let us now turn to the second assertion. Clearly, Rs] is Hilbert-Schmidt and self-adjoint. 
In addition, 

(i7°,,-c)V2(i[Q,i?,])±±(iyO^^_c)i/2 

and 

(ij;^,, - c)i/2i[Q, Rs] = i(i/0„. - c)i/2Qi?,- - iRsiH^,, - c^^^Q G ©2. 
Hence, i[Q,-R5] £ Q and \i[Q,Rs]\Q < 2\Q\q. We deduce from the fact that (see [3 Lemma 7]) 
VI < 39 < 00, yAe&p, limpA-A-Alle = 0. 

(54,0 

Let us finally prove the third assertion. We focus on the first estimate; the other ones can be 
established in a very similar manner. Consider for instance 

(1 - ^l,)V,RsH^,^Q{l - 70,,) = (y,i?5H0„0)++ = v;+RsH^,^Q++ + V^- RsH^,,Q'+ , 
the term -fp^j.VgRsHp^j.Q-fp^j. being treated similarly. Since Vg and RsH^^^. are bounded, Q^^ £ &i 
and i ^g"^ ^ ®2, the operator (Vgi?5iJpgj,Q)^^ is trace-class on L^(R'^). Besides, 

Tr ((y,i?,-i/0„O)++) = Tr(l//+i?5ifO„0++) + Tr(y+-i?,-i7°,,Q-+) . (73) 

The second term in ([75]) is the trace of Vg RsH^^^^Q^^ = Vg AsQ ^, where As = "f'^^^RsH^^^ is 
uniformly bounded in 5. It can therefore be bounded by C|| Vg+^ II 62 II Q II Si hence by C|| e|| i2p,c || Q|| Q 
in view of Lemma [J] 

The first term on the right-hand side of (|73p can be rewritten as 

Tr(v++RsHl^.Q++) = Tr(^*I,(i7°,, - c) V2g++(i/0^^ _ c) 



1/2 

where Ag = H^^^{H^^^ — c)^^Rs is uniformly bounded in 5 and w is defined in ([70)) . The bound- 
edness of w and the inequality ||i?5|| < 1 imply the existence of a constant c > 0, independent 
of 5, Q and Q, such that 

Tr ((y,i?,-i7°„Q)++) I < H||f?||L^ncllQlls. 
This therefore gives the expected estimate. □ 

5.4 Proof of Proposition [5] 

We start with the case n = 1. We easily deduce from Lemmas[51[3]and|l]tliat if p e L^(M+, L^(M^)n 
C) and if w = Vc{p), then Qi G C"(M+, Q) and the following estimate holds: 

VteM+, ||gi,„(i)lls</3(Ccom+/3Ccom.s||Q"||Q) / || p(s) || L^nC 

\p{s)\\L^-ncds. (74) 







< /3max (Ceom,/3Ccom,s) (1 + ||Q°||q) ^ 



We also infer from Lemmas [U [3] and S] that Ti'Q{Qi.v{t)) = for all t e K._|-. Still using those three 
lemmas, we obtain by an elementary induction argument that for all n > 2, Qn.v G C'*'(IR+7 Q)j 
Tro(Q„,t,(t)) = for ah < e K+, and 

VteR+, ||Q„.„(t)||Q</3Ceom.S / \\pis)\\L2nc\\Qn-l.vmQds. (75) 

Jo 

The estimate pop being true for n = 1 in view of (j74p . it remains true for all ri > 2. 

The right-hand side of (j27p therefore normally, hence uniformly, converges in Q on any compact 
subset of R+, to some Q{t) such that Q(-) G C°(K+, Q). It is then elementary to check that Q{-) 
is the unique solution to (E31)-(I11D in C°{R+, Q). 
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5.5 Proof of Proposition [6] 

We consider the regularized operator Xq based on p7p . and defined as 



V I— > 



The function f : q fq{-) is in L°° {T* , L'f^^j.{T)) as soon as 



„ +00 
sup f ( 



n,m— 1 



and 



|L°°(r*,L2^,(r)) < I sup -r ^ i'i-n<N<m + 'i-m<N<n) \0im,n,q,q'\^ dq' \ 



1/2 



(76) 



We show in this section that this operator is in fact well-defined and bounded from C) to 

L'^{R, L2(j^3) p q-j foj. g^j^y ^ > 0, so that S"'i is a bounded operator on L'^{R, L'^{R^)). 
In the sequel, we will meet expressions of the form 

„ +00 

fqi^)^'f y2 i'^n<N<m ~ ^niKNKn) am.n,q,q' Um.q'{x)Un.q+q'{x) dq' . 



(77) 



It is easily checked that the coefficients am,n,q,q' in the expressions below satisfy ([77]) using the 
following estimates. 



Lemma 19. 

(1) There exists (a_, 6_) e M:^- ^ '^'^'^ (0+7^+) e IR+ x K such that, for all q eT* , 

a^n^^^ + b- < £n,q < a+n^^^ + b+. (78) 

(2) There exists a constant C G K.-|_ such that, for any function v G H'^{R?) and for all 1 < n < N , 
m> N + 1, q eV*, 



\^n,qi ^q—q' ^m,q' / 



dq' < C\\v\\jj2m 



2 -4/3 



(79) 



(3) There exists a constant C G such that, for any K G 



(u„,„e'^-"M™,,.)i2 < C(l + |A-|^)m 



-2/3 



(80) 



Proof The bound ^ follows from dUD (see also (3.9) in [5]) and the resuhs of Section XIII.15 
in [TS]. To prove ((7^ and ((5(1)) . we rewrite, for m large enough, Um,q' as £^g>iHpc^)q'Um,q' , so 
that, for aU w G i/p(,r(r). 



1 



((-^^pcr)g'(«^'«n,9),Um,9') 



1/1— — — /I— — 
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We infer from ([75]) that there exists a constant C E M+ such that for aU 1 < n < TV, to > TV + 1, 

q e r*, q' e r*, w G h^,A^), 

I {Un,g,WU„i,q')^2 I < C|1w|1h2 m"^/^ 

per i'^' 

Choosing w(a;) = leads to ((80)) . Applying the square of the above inequality to u; = "g-g' 

for V e i/^(R''), and integrating on F*, we obtain 



(u„,„V,_,,U™,,,)i2^^ dq' <c(f dg' ) 771-4/3 < C'Ml^ TO-4/3 



which completes the proof of the lemma. 



□ 



The proof of Proposition |5] is performed in two steps: (i) we first give the expression of 
[J-'t{xov){uj)]q(x) since this quantity is the basis for several computations in this section and the 
following ones; (ii) we then evaluate (/2, (o^ fi) i^2(^]^2y The proofs are written for regular functions 
V, /i, /2, the general result following by the continuity of on L^(R, L^' 



Lemma 20. For any function v e ,y{R x M'^), the following equality holds in L°°{R x F*, ip(,j.(r)) 
(with w e M and q £ T*): 



^ ^ (ln<A^<m ~ ^'rn<N<n 



)— Z Z -^^^^ dq'. 



(81) 



£n,q+q' — £m,q' — — ITj 

Proof. Let w e ^(M X R3)_ We first note that 

Ql„{t) = -i (Uoit - sHsh'pMt - s)* - Uo{t ~ s)7°,,7;(s)f/o(< - s)*)e-''(*-^) ds 

[Uo{t - s)7°„.«(s)(7p°„.)^f/o(t - s)* - Uo{t - s){^"^,,)^v{sh"^Mt " s)*)e-''(*-^) ds. 
The Bloch decomposition of the operator Q ^ (i) reads 



9.9 



t/o(t-s)7°,, 7;(s) (7°,,)^[/o(i-,s)' 

J g L J q-q' L 



7°ert/0(t-s)* e-^^'-^^ds 



, L J g L J 5 

/. TOO 

= i / X! (ln<JV<jn - lm<Ar<Ti) - s) /l„,5,m,<j'(s) |u„,g)(-U„i,5'| ds, 

•'^ n,m=l 

where = exp (7] + i(£„,, -e,„,,'))i^ 1(>0 and hn,q..m,q'{t) = {Un^q, [v(t)]q^q>Urn,q') l^^^^- 

It can be checked that [Qi = /p. [*9i t, (0] ^ ^/ '^9' is well defined in Lp^j.(r) when r; > 

since gnqmq'i^) Uniformly integrable in t,n,m,q,q' . Therefore (see [SI Section 6.5]), the 
following equality holds in L°°{M. x F*,ip„(F)) for the function {t,q) i-> [xov{t)]q: 

x'Mt)] = / iQUt)] ,{^.^)dq' 

J 9 Jr* 9+9 ,9 

+ 00 

^ ^ (ln<Af<m ~ lm<Ar<ri) (^ffn.g+ij' ,m.g' ^n.q+q' ,m,q'^ {t) Ujn^q' {x)Un.q+q' {x) dq' . 
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Remark that ^, and hn^q+q\m,q' are both mtegrable, and that 

{9lq + q',m.q') i^) = ^ ■ 

It foUows that 

L J q 

+00 



7 X] (l"<Af<™ ~ lm<Af<n)-7^t (5r!,g+q',m,g') ('^) -^t i^n.q+q' ,m,q') (uj) Um.q' ix)Un^q+q' (x) dq' 

\ ^n,q+q' 1 (^) 



+00 

^ ^ (ln<7V<m ^ lm<7V<n)" 
n.m— 1 



■dq', 



where the equaUty holds in Lp^j.(T) (for functions in the x variable) uniformly in a; € M and 



ge r* 



□ 



Lemma 21. For any /i, /2 e ^{R x M^) anrf ?? > 0, 

(/2,(S"'/i)l2(L2) 

Tr 



-— Im 

TT 



(82) 



Proof. Using Lemma [501 

(/2, <§"'/l)L2(i2) 
+ 00 

Ed 



r* 



^t«y'(/2)(c^) (x) J-a^y'(/i)(c^) {x)dxdqdio 



+00 ^U,ri,g', [-^tWy^(/2)(w)]gWn 
n<N<m ^ '^m<N<n) 



g '^7l,q+q' 



Un,q+q', [Ffll^J'^ {fl){Lo)]q U„i.q' 



(r 



n.m—1 



-— Im Trr2 

TT \ P° 



(7° )^ 



(r.)2 (i?o„ + w + i?/)g' 



[-^t«y'(/2)(a;)]. 



£n,9+g' — £m,q' — W — 1?7 
(7pcr)g+9' 



■ dg' dq du 



[Ftvl'^{h){uj)]qdq' dqdiudz 



where we have used the fact that the terms in the sum over 1 < m < < n are the complex 
conjugates of the corresponding terms in the sum over 1 < n < N < m. Remarking that 



= TrL2 



(7° 



pcrJq' 



IT' Jt' z- {HO^^. + oj + iri)q, 



per 
(Tpcr) 



(7p°c..), 
z - {HO 



z - (H^cr + ^ + iv) 



— [Ttvl/'{.f2){-io)] 



z - H° 



por;<3 



■Ttv\l\h){^)\q-q' dq'dq 
dq' 



J q',q' 



Tr 



(7" 

V /pc 



' Z -"per ^ 



we obtain the expected result. 

Proposition ini now easily follows from ([5^ . using the density of x R^) in L^(R,i^ 
The bounds on (o^^uj) are a consequence of (see the discussion after this equation). 



□ 
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5.6 Proof of Proposition [7] 

Note first that, thanks to Lemma [T51 above, the smus over n, m in (|45p are convergent when 77 > 0. 
In addition, for all 77 > and all g <E F*, the expression ps)) can be rewritten as 



-iK 



X ('7por)g+9' AK'x 



■ dq' dz 



+ 



Z - {H^ei-)q+q' ^ - (i/0„ + UJ + 1??)^ 



Z - (-ffpor ^ - (^fpcr)q 



(83) 

where "^jj is plotted in Figure [TJ This gives the continuity of the mapping {lo, q) i— > x'i^^ l) "-"^ 
M X F* for any ry > 0. 

To prove (|31]), we use Lemma BOl and write 

{2i,f/^FtAx>){L0,q + K)= f \Wov){^)\ {x)e-'^-^dx 



+00 {^it„,,+5', [^"(^(w)]^?/^,,.^^^ ( 

= ^ ^ (lTi<Af<Tn ~ l?ri<JV<n) 



n,m— 1 



^771, q' 7^ ' ^n-q+q'"} 



Now, 



^ 71, q-\-q' ^m,q' ^ 1^ 



■ dq' . 



which implies that 



J^tAXov)iu;,q+K)= J2 



K'eTl* L n,m=l 



with 



£n,q+q' — £m,q' — W — 1?/ 



■dg'. (84) 



Therefore, 



+ 00 



Since F{vV^ f){k) = -v/lTrl/cl ^Tf{k), we obtain that the entries of the Bloch matrix of the 
operator S'^ = ^y^Xo'^c^^ are given by dH]). 

5.7 Proof of Proposition [5] 

The outline of the proof is the following: 

(i) we first characterize the limit when rj goes to zero of the matrices S'^ ^, for a given pair 
iK,K') e (7^*)2 (Lemma 1221); 

(ii) next, we show that for any / G S^{M. x M'^), the series 



'K.K' 



K,K'en* 



converges to J^t,x{S'^ f) in 5^'{M. x R"^), and has a well-defined limit when 77 goes to zero 
(Lemma [221) ; 
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(iii) finally, we prove that (^"^ strongly converges to S on time intervals of the form (— oo,r], 
which allows us to identify J^t,x{'S' f) with the limiting series obtained in the previous step 
(see Lemma and the discussion after its proof). 

Lemma 22. For any K,K' G TZ* , the family of functions S'^ x' defined by (|44p has a limit in 
X M.^), denoted by Sk,k' , when rj goes to zero. Moreover, the support of <§k,K' is contained 
inRxT*. 



Proof. It is easily seen that for any ?7 > 0, the function (w, q) i-> ^, (w, q) belongs to i°°(Rx M'^), 
hence to ^'(K x K^), and that its support is included in M x V* . Fix a function Lp e y{R x M^). 
It holds: 



/pii 





1 

W\ 



+ 00 



<N<7n ~ ^m<N< 



n) 



n,m— 1 

where T^^^ ^ ^/ is defined in ([84|l . Now, 



\q + K'\ 
\q + K\ 



'^Irn.K.K' , qM^ , q) duJ ^ {Um^g' , C 'Un^g+g 



-iK x , 



with 



pcr\ 



du). 



£n,g+g' ^ £m,g' — UJ — IT] 

Standard computations show that the functions v]/^.™''; are bounded in L°°{T* x F*) uniformly 
for < 77 < 1, and 1 < 71 < iV < m or 1 < to < iV < ?i. In addition, 



\im%-^'^'^{q,q') = ^"^'^iq,q'), 



where 



$"'"(9,9') :=-(p.v 



,Lp{en,g+g' - era,g' + ■,q)) - mipiSn^g+g' - Ern.g' , q) ■ (85) 

.y"..y 



It then follows from Lemma flQl that \<S'x k'^V) has a limit {Sk,k' i'^) y when rj goes to 
zero, given by 

^ +00 



Af<m ~ -'-m<Af<nJ 



\q + K'\ 
\q + K\ 



iKx \ / 

^n,g+g' / ji^2 \^n,g+g' : 



We also infer from the above arguments that there exists a constant C independent of Lp, K and 
K' such that 

{<^K,K' .f^^) y, y \<C{l + \K'\f{l + \K\)^r'{^l 
where the seminorm .yVr* is defined on 5^(R x M^) by 



^r'{^) :=||(l + |c.|)^|| 



L=^(Rxr*) 



dip 



The limit S'k.k' of K' therefore defines a tempered distribution of order 1. Besides, as the 
distributions (S*^ ^, are all supported in R x F*, so is their limit. □ 
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Lemma 23. Let f e ,5^{M. x M?). For all tj > 0, it holds 



(86) 



in X M^). In addition, the above quantity converges in x M'^), when rj goes to zero, to 

the tempered distribution 



where the tempered distributions <Sk,K' o,re defined in Lemma[ 



Proof. The computations performed in the proof of Lemma [52] show that there exists a constant 
C > and 770 > small enough such that, for all < ^/ < 'yo, all K,K' in TZ* , and all ip G 

X M.^), 



<G{\^\K'\Y{\^\K\)^^,{^). 



Therefore, 



< C{1 + \K'\f{l + \K\) ^r'{T.K'{J't,j) t-k{^)) 

< C{1 + \K'\fil + \K\) J^r- {{J't..f){; ■ + K')) ^r- 

Consequently, for all / € o5^(]R x R'^), and any < 77 < 770, the series 



converges in yiM. x R^) and 



where for (p, q) G N x N, ^p^q denotes the Schwartz seminorm on 



defined as 



max 



ii + \io\r{i + \k\y 



\o!p\<P, \aq\<q 

The claimed convergence result is then easily obtained. 



□ 



Remark 24 (Sufficient regularity requirements on the function /). The above proof shows that 
the series (|86p are well defined as soon as 7(^-4. x/) < 00. Actually, weaker conditions such as 



\\il + H)il + \k\)^+'Tt,J{uj,k)U^ + 
can be derived by using sharper estimates in the above two lemmas. 



< 00 
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Lemma 25. For all f e L'^ {R, L'^ (R^)) , it holds 



VTeM, lim(f'' = (f/ in L^{{-(X),T],H\R^)). 



Proof. This result is a straightforward consequence of the foUowing fact: for any given potential 

V e L\R,C'), 

VTeM, sup \\Ql,{t)~Qi,M\o (87) 

te{-oc.,T] ■ iiy ,,->o 

together with the continuity of the hnear mappings Q 3 Q ^ pq e L^(M^)nC, vl^^ : L^(]R^) — > C 
and t-y^ : L^j-j^S) p ^ ^/ p L2(E3) = H^{R^). Actually, it is sufficient to show ^ for 

V G C^(R X M'^). Indeed, fix e > 0, and approximate v by some v G C^{R x R^) in such a way 
that 

\\^it)-v{t)\\c dt<s. 

R 



Then, using Lemmas [5] and SI we obtain 



ft 

Vte 



and similarly \\Qi,vit) — Qi,vit)\\Q < (3Ccom£- Let us now consider v € C^(M x R^) and To > 
large enough so that the support of 1 1— > v{t) is contained in [—Tq, Tq]. Then, 

sup ||Q^(i)-Ol,.W|L </5^com|k||Li(K.C') (e"^'^'+'^''^-l) — >0, 

which concludes the proof. □ 

With these results, the proof of Proposition |S] is now straightforward. Indeed, Lemma [55] 
implies that, for any / € y{R x R^), Tt,x (^''/) converges to Tt,x (^/) in ^'(R x R^), while 
Lemma [221 allows to identify the corresponding limit. 



5.8 Proof of Proposition (TU] 

We fix / g ^(R X R3) and prove that Tt,x (^"./) converges in ^'(R x R^) to Tt,x (^°/) when 
a goes to zero. 
The expression 

QIM = { uo (^) WUMs)] uo ds 

shows that the adiabatic evolution can be understood as the standard evolution upon considering 
the evolution operator with generator H^^^/a (hence replacing en,q+q' —£m,q' by {Sn,q+q>—£m,q')/ct 
in the expressions involving Bloch matrices), and rescaling globally the result by a factor a~^. 

According to Proposition [5| and the results established in Section [5771 the quantity Tt^x (^S""/^ can 
therefore be expressed in terms of the Bloch matrices {(a^ Ki)K,K'ew by the following equality in 

y"{R X r3). 

where for any ip € y{RxR^), 

+ CXD 

('^K.ifM t/f^^^ ^ = ^ (lri<JV<m - lTn<Af<n)£']*.'']^'"(V?), 
' n,m— 1 
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with 



and 



ir"fa»'):=4(p.v.(^),.( -'--'^'7"--'^ H--.,)) 



, ,q 

On the other hand, 

where for any ip € .y{R x R^), 

+00 

' n,m— 1 

with 

We now use the same arguments as in the previous section to prove the convergence of k'i'^^ 
to (^Si( j^i, (p^ when a goes to zero. The only point we need to check is that 

• there exists a constant C such that for ah < a < 1, all 1 < n < iV < ?7i or 1 < to < < n, 
and all q, q' in F*, 



• for all 1 < n < < m or 1 < 7Ti < < n, and all q, q' in F*, 

lim$;j'™^"(g,q') = ^ / ^{u:,q)duj. 

"4-0 en,q+q' ^ £m,q' JR 

This is a direct consequence of the following lemma. We denote by y^p, p E N, the usual Schwartz 
seminorms on S^{M.). 

Lemma 26. For given ip £ ^(M), ?/ G M, and a > 0, we set 



Then, there exists < 00 such that for all \y\ > 6 > 0, and all a G (0, 1], 
and, for any y ^ 0, 
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Proof. Consider the case when y > 0. It holds 

h4,,y{a) = ^g^ (I) , 



where 



9'4'iz) 



dx ^ z 



^ ij){z + x) ~ ijj{z — a;) 



dx — z 



+00 



The third term can be bounded as follows: for z > 1. 



+00 



-0(2 + X) 



dx 



z {z + x)'^%Ij{z + x) 
z + X x{x + z) 



dx 



< 



+00 



dx 



x{x + z) 



ip{z + x) 



dx. 



0. 



Moreover, for 2; > 2, 

'iIj{z + x) — ip{z — x) 



dx 



= 2 



zip' {z + Ox^zx) dx 



< 



2z 



where O^.z G [^1; !]• Lastly, for z > 2. 



tlj{z + x) 



dx 



1 



-%j}{y) dy 



(z-\)l2 



1 



1 - yjz 



{z-\f 



ip{y)dy- 



z^-\-oo 



^ 0, 



1 - y/z 

The first term of the right-hand side of the above equality is controlled by 

'■(z~l)/2 



{7.-1)12 ^-y/z 



ipiy) dy. 



IV'I <2 



1 + 



and converges to J^ip when z goes to +00. To bound the other term, we notice that 



{z-i)/2 ^-y/z 



■0(y) dy 



< z 



\i^iy)\dy < z 



{z-l)/2 



{z-i)/2 r 



■ dy - 



2(z-l)2 



2:— ^-f Ctt 



^ 0. 



Hence the function z i-> gip{z) is bounded by C^3{'ip), uniformly on [2, +00), and converges to 
/jgV' when z goes to +00. This completes the proof. □ 

5.9 Proof of Theorem 1111 

To simplify the notation, we denote by v{s) := —Vc{i^{s)) and w'^{s) := v{s) + Vc{pq(s)) = 
Vc{pq(s) — i^{s)). We proceed in three steps: 

(i) we first show that the dynamics is well defined for short times fSection lS.y.ip : 

(ii) we then extend the result to arbitrary times using some energy estimate fSection l5.9.2p : 

(iii) we finally establish a few qualitative properties of the solution fSection l5.9.3p . 

5.9.1 Local-in-time existence and uniqueness 

The existence and uniqueness in C"([0,T],Q) of the solution to the integral equation ([55]) for 
short times easily follows from Lemmas |2] and SI using standard techniques. For each T > 0, 
we consider the mapping Ft : C°([0,T], Q) C"([0,T], Q) defined by 

FTiQm^Uom^'Uoitr ~i f Uo{t~s)\v,{pQ^,)~iy{s)),^^^^ + Q{s)\Uo(t-srds. (88) 

Jo L J 

Notice that the solutions of the integral equation ([55)) on [0,T] are the fixed points of Ft. For 
each T > 0, the mapping Ft is well defined in view of Lemmas [H [2]and|ll The existence of a 
fixed point for T small enough is, in turn, given by the following 
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Lemma 27. For any R > I3\\Q^\\q (where j3 is defined in Lemma[^, there exists T > small 
enough such that Ft is a contraction on 



sup \\QmQ<R 

0<t<T 



Proof. Let Q E Bf>. Lemmas [U [3] and |H together with the continuity property show that 
\\FT{Q)it)-Uom°Uo{tr\\Q<p(c,omJ^ \\pQis)'iy{s)\\cds + C,on..Ql^ \\pQis)-iyis)\\L-nc\\Q{s)\\Qds 

< /3(Ccom + C'com.Qfl) / \\ PQ{s) - l^is)\\ L^nC ds 

is)\\L^ncds . 







Therefore, 

WFriQmWQ < /3||Q"'||s + /3(Cco,„ + Ceo„.,Qi?) I CpRT + I Ms)\\L^.ncds I < R 



cT 







for T small enough, so that the application Ft : Bp — Bn is well defined. Now, consider 
Qi,Q2 G Br. Then, 



FT{Qi){t) - FT{Q2m = -\ I U„{t-s) [u;«i(s)-u;«^(s),7°„, + gi(s)] [/o(t-s)*ds 

Uoit-s) [u;'3^(s),Qi(s)-Q2(s)] Uoit-s)*ds, 







so that 



rt 

iFriQim - ^^T(Q2)(t)||Q < /3Cp(Ceom + Ccom^Qi?) / ||gi(.s) - Q2{s)\\Qds 





+ /3Ceom.Q / i\W{s)\\L-^nC + CpR)\\Ql{s) - Q2{s)\\Qds 
JO 

<a(T) max ||gi(s)-Q2(s)||Q, 

u<s< i 



where 

r-T 



a(T) = /?Ceom,Q / |k(s)||L^nC + /^C-p (Ceom + 2Ceom,Qi?) T 



is (strictly) smaller than 1 when T is small enough. This shows that Ft is a contraction provided 
T is small enough. □ 

5.9.2 Global-in-time existence and uniqueness 

Let [0, Tc), < Tc < oo, the maximal interval on which the solution to the integral equation (|53p 
is well-defined. In order to obtain global-in-time existence and uniqueness, that is to prove that 
Tc = oo, it suffices to show that ||Q(i)|lg does not blow up in finite time. For this purpose, we 
rely on the following energy estimate: 

Vte[0,r,), £(t,Q(t))=5(0,g")- /" Z?(pQ(,),i/'(s))ds, (89) 
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where we recall that £{t,Q) is defined m ([M)) as 



Sit, Q) ^ Tro(i/p",,g) - DipQ,i.{t)) + -D{pq,pq). 



Although the formal derivation of (|89| is a simple exercice, the rigorous proof is somewhat tech- 
nical. We first complete the proof of the global-in-time existence and uniqueness, assuming that 
([5n|) holds true; the latter equality will be established at the end of the present section. 

From ([M)) . we infer that 

£{t, Qit)) < £{0, Q°) + 1^ |i?(pQ(,,), i^'(s)) I ds < £{0, QO) + Cp \\iy'is)\\c\\Q{s)\\ q ds. 

On the other hand, we deduce from Corollary 2 in [S] that there exist a,b > such that for all 
t G IR+ and all Q e /C, 



£{t, Q) > allQIls -I- eFTro(Q) - b - -||Ki)llc- 



In addition, Tro(Q(t)) = Tro(Q°) for all t G [0,Tc) in view of Lemmas [H [3] and S] and the formula 



Q{t) = Uo{t)Q"Uo{tr - i / Uo{t - s) Vc{pQis) - Hs)),Jpor + Q{s) Uo{t - sYds. 



Therefore, 



where 



Vie[0,T,), ||Q(i)||Q<ai(t)+ / «2(s)||Q(s)||sds 



ai(i) = - (0,Q°) -eFTro(g") +5 + , = ^Hz^'WIlc 

a \ 2 J a 

As, by assumption, v £ W^^^{M.+ ,C), we infer from the Gronwall lemma that does not 

blow up in finite time, which implies that Tc = oo. 

Let us finally establish the energy estimate (|5^ . The proof is based on the following result. 



Lemma 28. Consider the regularization operators Rs — (l + (5 |iJp(,j. — ep |) ; o,nd the regu- 
larized energy £s{t,Q) = £(t, RgQRs), where £ is defined in (j54p . Let v G WjqJ(R+,C) and 
Q G C^{[0,Tc), Q) be the unique solution to the integral equation i53\} . Then, there exists a 
constant c G M+ such that, for all t G [0,Tc), 



£s{t,Q{t)) = £s{Q,Q'')- [ D(^pR^Q(,)R^,v'{s))ds+ [ rs{s)ds, 



(90) 



with 



\rsit)\<c[\WmL^^nc + \\Qm 



(1 + \\QmQ)\\RsQ{t)Rs - QWIIq + mQit),Rs] 



RsQ{t)Rs - Qit) 



(91) 



where Q{t)=i[wQit),^°p^^. + Qit)] 



The energy estimate (|89p can then be easily deduced from (j90p , by remarking that the mapping 
Q 3 Q 1-^ £{t,Q) E M. is continuous, and that the first two assertions in Lemma \TE\ allow to pass 
to the limit in (1901) by means of the dominated convergence theorem. 
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Proof of Lemma \28l In this proof, the constant C > may vary from hne to hne, and can be 
chosen to be independent of 6. By density, it is enough to estabhsh ((90)) in the case when ly S 
C\R+,L%R^)r\C). 

The solution of ([5^ can be rewritten as 

Q{t)^Uo{t)(^Q" -ij^ Uoisr[w'^{s),^^^^ + Q{s)\Uois)ds^Uo{tr. (92) 

In addition, 

where RgHp^^. — H^^^.Rs is bounded. Therefore, 1 1-^ RsQ{t)Rs is difFerentiable and 

ij^(^RsQ{t)Rs) = RsH°^Mt)Rs - RsQit)RsH^,, + Rsiw'^ (t) n^cr + Q{t)]Rs- (93) 

Since Q{t) G Q and \[w^{t),^^^^ + Q{t)\ G Q by Lemmas [3] and SI it is easily verified that 
ft{RsQ{t)Rs) e Q for aU S > 0. Now, ^ implies 

TYo{H^,,RsQ{t)R5) = Tyo{H^,,RsQ°Rs) - f Tro {H^,,Rsi [w'^{s), 7per + Qi^)] Rs) ds, 

Jo 

so that 1 1-> Ti-o{Hp^^.RsQ{t)Rs) is differentiable, with 

j^(Tio{H°,,RsQ{t)Rs)) = -Tro {H^.rRsi K(t),7pcr + Q{t)] Rs) 

= -Tro {H^^rRsi [wQ{t),Q{t)] Rs) 
= -Tro {RsH^^^Rsi [w^ {t),Qit)]) , 

where we have used that '-fp^^, Rs and RsH^^^. commute. Moreover, 

^^D{pR^Q(t)R„v{t)) ^ D[pji^Q(t)R,,v'{t)^ +Tro (^^{v{t))^^{RsQ{t)Rs 
where the second expression makes sense since ^{RsQ{t)Rs) G Q. Finally, 

^^D{pF.^Q(t)R,: PR,Q(t)Rs) = 2Tro {vc{pR,Q(t)Rs)^^{RsQ{t)Rs) 
Therefore, t ^ £s{t, Q{t)) is differentiable and 

j^(^£sit.Q{t))) ^-D(^pR^Q(t^R^,v'{t)) +rs{t) 

where 

rs{t) - Tro [w''''^''' {t)j^(RsQ{t)Rs)^ - Tro {RsH°^^Rsi [w'^ {t) , Q {t)]) . (94) 
Therefore, ([50]) holds true for rs given by ([M]). Using ([M)) . we may rewrite rs{t) as 
ra(i) = -iTro(wc(pfl,Q(t)i?,-g(t))i?5-ffper<3(i)^<5) + T^T^T^o[MPRsQit)Rs-Qit))R5Qit)RsH^cr 



Tr 



(^;c(Pi^,Q(t)i^,-Q(t))i?^i[^«^(^),7pcr + Q(^)]^^)) - Tro (t)i?5i[«;^ (t), 7per + 



- Tro {wQit)RsH°^^i[Q{t), Rs]) - Tro (w'^mQit), Rs]RsH^^^) . (95) 
We deduce from the third assertion in Lemma [18] that 

TroUipR,Q^t)R,_Q^t))RsH'p,,Qit)Rs) < C \\Q{t)\\Q\\RsQ{t)Rs - QmQ, (96) 
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and 

TTo[v,{pR,Q^t)R,^Qit))R5Q{t)RsH°,,^\<C\\QmQ\\RsQit)R5~Qm^ (97) 
It also follows from Lemmas |3] and |4] that 

Tro (wc(pfl,Q(t)fl^,-Q(t))-R'5i["''^W'7pcr + Q{t)]R5^ 

< c (1 + IIQ(t)lls) (Ik(t)llL^nc + IIOWIIq) WRsQims - QWIIq- (98) 

To bound the fourth term in the RHS of ([55)) . we notice that 

Tro {w'i{t)RsAw'i{t),^l, + Qit)]Rs) - Tro («;^(t)i[«;Q(0, 7p°c.. + Qit)]) 

+ Tro {w^{t) [RsQ{t)Rs - Q(t))) , 

where we recall that Q{t) := i[w'^ (t) , + Q{t)]. Now, for ah V € Vc{L'^ n C) and all Q € Q, 

Tro (l^i[y,7pcr + Q]) =0, 

so that 

|Tro(«;Q(i)i?5i[i«'^(0,7per + Q(i)]fi5)| <C (|li.(0llL2nc + |lQ(t)||Q) \\R5Q{t)Rs - Q{t)\\^. 

(99) 

We finally infer from the third assertion in Lemma [TS] that 

{Troiw'^msH^eAQW.Rsdl < C (11^^(^)1^^00 + 11(3(0112) UQ{t), RsWIq, (100) 

and 

\TTo{w'imQ{t),Rs]RsH^,,)\<C {\\u{t)\\L2^c + \\QmQ) MQ{t), Rs]\\q. (101) 
Collecting (|96| - (fT0T|) . we obtain (|9l|). □ 

5.9.3 Properties of the solution 

The preservation of the trace has already been proved at the beginning of Section [5.9.21 

Let us now assume that 7°^^ + is an orthogonal projector. Since 7°^^ + Q{t) is self-adjoint 
and non- negative, proving that 7pgj. + Qit) is an orthogonal projector amounts to proving that 
(7pcr + Q(i))' = lpe. + Qit) for all t > 0. Introducing r(i) = h^,, + Qit)f and 7(f) = j^,, + Qit), 
it formally holds 

i|7(t) = [i^(t),7(i)], 
where H{t) = iJ^^, + Wc(/5Q(i) - !^(i)), and 

ijTit) = [i7(t),7(i)]7(i) +7(t)[i^(0,7(i)] - [i^(t),7(0'] = [^(0,r(t)]. 

The above formal computation can be made rigorous upon using mild formulations, and estab- 
lishing the second equality by some limiting procedure involving regularization operators, as in 
Section [5.9.2l We do not detail this point here for the sake of brevity. The uniqueness of the mild 
solution of the linear equation 

and the fact that 7(0) = r(0) allow to conclude that 7(t) = T{t) for alH > 0. 
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5.10 Proof of Proposition [T3l 

Proposition [5] shows that the regularized operator 

(with Xo defined in ([75]) ) is such that, for any q\,Q2 G 

(£»2,/^''6Ii)l2(C) 

(7p°c..)^ 



-Im 



Tr 



per 



\ri) 



(7p°o..) 

^ -'-'per 



dijj dz 



(102) 



(103) 



where the contour is plotted on Figure [T] 

We now investigate the limit 77 — > of the latter expression. To this end, we choose a contour 
similar to the one of Figure[2j such that, for all i] E [—1, 1], 



dist(^, (a n [£f, +00)) ± O + ir;) > > 0, 



(104) 



and consider gi,g2 <E x M.^) such that J^t0i,J^tQ2 have support in [— r2,r2] x R^. We can 

then pass to the limit 77 1 in (|103p . and the limit actually make sense when TtVc{gk) G i^(R, C) 
for fc = 1,2 and J^tVc{.Qi)^tVc(,Q2) = outside a compact subset of {~g,g). This is the case 
in particular when both functions gk belong to Jifn (note however that it is even possible to 
give a meaning to this expression when only one of the functions is in the other one being 
in L^(]R,C)). The resulting expression is clearly symmetric in gi, g2'. 



Tr 



\ /per/ 



Z-(FpV + ^) 



TtVc{g2){i^) 



(7p°c..) 

^ ^^pcr 



TtVc{gi){uj) 



duj dz 



This finally shows that C, restricted to C L^(R, C), is a well-defined symmetric operator. It 
is also clearly bounded, with a bound proportional to {g — Jl)^^ in view of (|104p . 
Besides, for g € Jfa, 



l<?i<Af<?ri 



n 



{Un^q+q', [ J"t Wc (fi*) (w)] , U. 



m.q' / 1^2 



■ dq' dq duj > 0, 



-n,q+q' 



which shows that C is nonnegative on Jifii when < < g. 

Remark 29. Note that the corresponding expressions of {Cgi,g2)c ^.i^d {Cg,g)c in the time- 
independent (static) setting (see the proof of Proposition 2 in J^) are recovered by removing the 
integration over lo G [— f^, fi] and setting a; = everywhere. 

5.11 Proof of Proposition 1141 

This proof strongly relies on the proof of Theorem 3 in [6] and we only present the required 
modifications. For the ease of comparison, we follow the notation of [6] and hence use 77 for the 
spatial dilation operator 

{Ur,h){t,x)^T^^/^h[t,T^x). 

Computations similar to the ones performed in [SJ Section 6.10] give 



u;A~^u^v^f\ 

where A = ^^^(l + C)vc is bounded and invertiblc on the space 

^2 = {/ e i'(K, l2(m3)) I supp(j-j) c [-n, n] X 



(105) 
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^ ^ \/'2 1/2 

and where v and J"^ are the functions of Mti defined as = Uc v and P = Vc ' W^, with 
= ^U*Vc {yn — Pv^)- Note that Urj is unitary on H^i, with adjoint {U,j)* = U,j-i. The proof 
of the resuh therefore amounts to identifying the weak hmit of U*A^^U,^v. 
Lemma [20] shows that, for any function h G -L^(]R, ip(,j.(r)), 

Ft{Aqh){uj,x) = h{ui,x) 

,1/2 



- > (ln<W<m - lm<Ar<n) 7 (^'c ) q[Un,q+q' U,n,q']ix) dq . 

, JT' £n,q+q' ~ £m,q' ~ W 

This defines an operator A{uj, q) acting as (Aqh) (w, •) = A{uj, q)Tth{io, •), with (for g G Lp^^{r)) 

Aic.,q)g = g + {vl/')J ( (.y^),g W^g^ 

"-^ ^ ^ (106) 

r ^ / h° )^ , \ 



As A~^ is bounded on the weak limit of the left-hand side of pOSp can be studied for func- 
tions whose space-time Fourier transforms have compact support. For hi and /12 in L^(M, L^(M'^)) 
with compact supports, and for 77 small enough, 



{u;A-^U^hi,h2)^,^^,^ = (j-t(A-iC/„/ii)^(c^,.),-Ft([/„/i2),(w,-))^^ dqdoj 

JA{LL>,r]k)]~^ eo,eo^^^ Tt^xhi{uj,k)Tt.xh2it^,k) duj dk, 



where we have used that, for h e L^(R, L^(M^)) 

(2^ 

|r|i/ 



(■Fffe)g(a;,x) = ,^ Tt,xHuj,q + K)eKix), 



so that 

j:,{U,hUu.,x) = r"'^^^ E ■^.../^(-,^)eK(x). 

I I Ken' ^ '/ / 

Note that, since the spatial Fourier transforms of the functions hi, h2 have compact supports, for 
7; small enough, only the first component K = remains in the above sum. 

It is therefore sufficient to understand the limit l\^A{uj, Tjk)] ^ eo, eo^ when 77 — > 0, with eo = 

|r|-i/2. This is given by the following lemma, whose proof is omitted since it is a straightforward 
modification of Lemma 6 of [6] (note that Eq. (69) in is replaced by (jl06l) '). 

Lemma 30. Denote by Po the orthogonal projection on {eo}^. The following hold: 

(i) For all uj G {—g,g) and all a G A(LU,ria)eQ converges strongly in Lpgj.(r) to h„{ui,x) as 
7; — > 0, where for all k G M'^, the periodic function &fc(a;, •) is defined by 

fefc(w,-) = (|fc|2 + fc^L(w)fc)eo 



(107) 
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with 



N +00 



£n.a - W 



and where Gq is the operator defined on L'^^j.{T) as 



KeTC*\{o} ' ' 
and which satisfies PqGq^'^ ~ Gq^'^Pq. 



\T\i 



iK-. 



where fx = f{x) -|^^ji77 dx, 



(ii) For a given lo G {—g,g), the family of operators PQA{Lj,q)Pf), seen as hounded operators 
acting on PqL^^^{T) , is continuous with respect to q and 



PoA{uj, q)Po 



strongly as q ^ 0, where G{lu) > 1 is the bounded operator on Poipcr(r) defined by 



C{uj)f = f+Gl'^ 



N 

i ^1 — 1 



(7" )^ 
V ipcvJq' 



(7° 



cv)q' 



(^pcr)?' ^n,q' ^ (^pcr)?' ^n.i 



for all f e Poiper(r). 
(Hi) For all cr G §^ and lo G {—g,g), the following limit holds: 



lim (eo, [A{uj, r]a)] ^ eo) = ^ 

^ l + <7TL{oj)a-(^Pob,{u,-),[C{io)r'PobA^,-)) 



Un,q' dq' 

(108) 



(109) 



We may now define the macroscopic dielectric permittivity as 

k^eM{i^)k = \kf + k^L{Lu)k - (^Pobkiu, •), [G{u)]-' Pobk{oj, •)) 



(110) 



The matrix inequality eM(w) > 1 is a straightforward consequence of p09D . using the fact that 
A>1. 
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